Low density ferromagnetism in biased bilayer graphene 
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We compute the phase diagram of a biased graphene bilayer. The existence of a ferromagnetic 
phase is discussed with respect both to carrier density and temperature. We find that the ferro- 
magnetic transition is first order, lowering the value of U relatively to the usual Stoner criterion. 
We show that in the ferromagnetic phase the two planes have unequal magnetization and that the 
electronic density is hole like in one plane and electron like in the other. 
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Introduction. — Graphene, a two-dimensional hexago- 
nal lattice of carbon atoms, has attracted considerable 
attention due to its unusual electronicproperties, char- 
acterized by massless Dirac fermions [l[ HJ ■ It was first 
produced by micromechanical cleavage of graphite and 
its hallmark is the half integer quantum Hall effect Q . 

In addition to graphene, few-layer graphene can also be 
produced. Of particular interest to us is bilayer graphene 
(BLG), where two carbon planes lay on top of each other 
according to AB-Bernal stacking. In BLG it is possible to 
have the two planes at different electrostatic potentials. 
As a consequence, a gap opens at the Dirac point and 
the low energy band acquires a Mexican hat dispersion 
!4|. This system is called a biased BLG, and provides 
the first semiconductor with a gap that can be tuned 
externally @, 0, 0| ■ Due to the Mexican hat dispersion 
the density of states (DOS) close to the gap diverges as 
the square root of the energy. The possibility of hav- 
ing an arbitrary large DOS at the Fermi energy poses 
the question whether this system can be unstable toward 
a ferromagnetic ground state - a question we want to 
address in this Letter. From the point of view of the 
exchange instability, BLG was found to be always unsta- 
ble toward a ferromagnetic ground state for low enough 
densities @,0|. 

The question of magnetism in carbon based systems 
has already a long history. Even before the discovery of 
graphene, graphite has attracted a broad interest due to 
the observation of anomalous properties, such as mag- 
netism and insulating behavior in the direction perpen- 
dicular to the planes 13, 11, 13, 3, 14 1. The research 
of s — p based magnetism (id. Il6l Il7| was especially 
motivated by the technological use of nanosized parti- 
cles of graphite, which show interesting features depend- 
ing on their shape, edges, and applied field of pressure 
[181 ] . Microscopic theoretical models of bulk carbon mag- 
netism include nitrogen-carbon compositions where fer- 
romagnetic ordering of spins could exist in ir delocalized 
systems due to a lone electron pair on a trivalent element 
[191 ] or intermediate graphite-diamond structures where 
the alternating sp 2 and sp 3 carbon atoms play the role 



of different valence elements [20|. More general models 
focus on the interplay between disorder and interaction 



2lLl22|. Further, midgap states due to zigzag edges play 
a predominant role in the formation of magnetic mo- 
ments I23l |24| which support flat-band ferromagnetism 




261 . |27| . Magnetism is also found in fullcrcnc based 
metal- free systems pjj ]. For a recent overview on metal- 
free carbon based magnetism see Ref. [29(. 

Model and mean field treatment. — Due to the electro- 
statically invoked band-gap, there is a large DOS fol- 
low carrier density and thus effective screening of the 
Coulomb interaction. Coulomb interaction shall thus be 
treated using a Hubbard on-site interaction. 

The Hamiltonian of a biased BLG Hubbard model is 
the sum of two pieces H = Htb + Hjj, where Htb is 
the tight-binding part and Hjj is the Coulomb on-site 
interaction part. The term Htb is a sum of four terms: 
the tight-binding Hamiltonian of each plane, the hopping 
term between planes, and the applied electrostatic bias. 

with 
r (r- 
and 



We therefore have Htb = S t =i Htb,,. + H± + Hy 
(r)[6«r(r) + bj(r - oi) + b L 



Htb,,, = -tJ2 r ,a a , 
a 2 )] + h.c, H ± = -t_i_Er J 



H v = lE r! x,,Kla(r) 



given by H v = U 
We used n XLa (r) 



a[ a {r)b 2a {r) + h.c], 
n X 2a(r)] . The term Hjj is 
"xiT { r ) n xii ( r ) + n x2\ (r )n x2 i (r)] . 



= -r-t 



xl a (r)x oa (r), x = a(b), as the 



number operator at position r and sublattice Ai {Bl) 
of layer i = 1,2, for spin a =1",!; ai = a(l,0) and 
a 2 = a(l, — \/3)/2 are the basis vectors and a ~ 2.46 A 
the lattice constant. Unless stated otherwise, we use 
t = 2.7 eV, t ± = 0.2i, and V = 0.05 eV [U. 

The problem defined by Htb + Hu cannot be solved 
exactly. We adopt a mean field approach, recently ap- 
plied to describe magnetic properties of graphene nanois- 
lands 3li]. Since the two planes of the BLG are at dif- 
ferent electrostatic potentials, we expect an asymmetry 
between layers for the charge density n and the magneti- 
zation m = n-f — n | (per unit cell). Accordingly, we pro- 
pose the following broken symmetry ground state, which 
also defines the mean field parameters: (n x i a (r)) = 
and (nxaa(r)} = + a 1 - 
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FIG. 1: (Color online) (a) Bilayer graphene DOS for U = 0. 
Inset: Zoom near the gap region, (b) U c vs Sn in the low 
doping regime. Inset: The same as a function of p. 



where An and Am represent the charge density and 
the spin polarization difference between the two lay- 
ers, respectively [32[. This leads to an effective bias 
V„ = V + C/An/4 - all Am/4. 

If one assumes the ferromagnetic transition to be sec- 
ond order, with m = and Am = at the transition, we 
are lead to a [/—critical U c given by, 



U c = l/pb(p,,U c ) ■ 



(1) 



where Pb(p, U c ) is the DOS per spin per lattice point and 
/i = ji — nU c /8, with chemical potential p. Although 
Eq. (JTJ looks like the usual Stoner criterion, the effective 
bias V a depends on U due to An. This makes Eq. (JTJ 
non-linear, and U c has to be found numerically in a self- 
consistent way. 

Simple results. — We start with the zero temperature 
(T = 0) phase diagram in the plane U vs Sn, where 
Sn is the doping relatively to the half filled case. An 
approximate analytic treatment is possible in this limit, 
which is used to check our numerical results. 

In Fig. [TJ (a) we represent the DOS of a biased BLG 
with U — 0. As seen in the inset, the DOS diverges at 
the edges of the gap. As a consequence, the closer the 
chemical potential to the gap edges, the lower the critical 
U c value. The low doping U c value - given by Eq. (JTJ in 
the limit UAn -C V - is shown in Fig. [TJ (b), both as 
a function of Sn and p (inset). The lowest represented 
value of U c is about U c — 2.7 eV to which corresponds 
Sn ~ 2.5 x 10~ 5 electrons per unit cell. The step like 
discontinuity shown in panel (b) for U c occurs when the 
Fermi energy equals V/2, signaling the top of the Mexi- 
can hat dispersion relation. 

It is clear from Fig. [TJ (b) that in the low doping limit 
U c is a linear function of Sn. To understand this behavior, 
first we note that for very low doping the DOS close to the 
gap edges behaves as pb(p) oc (\p\ — A s /2) -1 / 2 , where A g 
is the size of the gap. Using this approximate expression 
to compute the doping, Sn cx sign(/2) x /Jf' / 2 ^ x Pb( x )> we 
immediately get Sn oc sign(/i) / Pb{p) and thus U c oc |<5n|. 
In Fig. [TJ (b) both the numerical result of Eq. (JTJ and the 
approximated analytical result just derived are shown. 
The agreement is excellent. 
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FIG. 2: (Color online) Panels (a), (b), and (c) show the 
T — solution for m, Am, and An, respectively. Panel (d) 
shows the U vs Sn phase diagram at T — 0: symbols are in- 
ferred from panel (a) and signal a first-order transition; lines 
stand for the second-order one given by Eq. (JTJ. Labels: P- 
paramagnetic, F-ferromagnetic. 



Self- consistent solution. — In order to obtain the T = 
phase diagram of the biased BLG, we study how m, Am, 
and An depend on the interaction U, for given values of 
the electronic doping Sn. 

In Fig. [2] (a) it is shown how m depends on U for 
different values of Sn. The chosen values of Sn correspond 
to the chemical potential being located at the divergence 
of the low energy DOS, which explains the smaller critical 
U c value for smaller Sn. It is interesting to note that the 
saturation values of the magnetization correspond to full 
polarization of the doping charge density with m = Sn, 
also found within a one-band model i . In Fig. m (b) 
we plot the Am vs U . Interestingly, the value of Am 
vanishes at the same U c as m. For finite values of m 
we have Am > m, which means that the magnetization 
of the two layers is opposite and unequal. In Fig. [2] (c) 
we show An vs U. It is clear that \8n\ < \An\, which 
implies that the density of charge carriers is above the 
Dirac point in one plane and below it in the other plane. 
This means that the charge carriers are electron like in 
one plane and hole like in the other. As U is increased 
An is suppressed in order to reduce the system Coulomb 
energy. 

In Fig. [2] (d) we show the T — phase diagram in the 
U vs Sn plane. Here we concentrate on the V = 0.05 eV 
case. Symbols are inferred from the magnetization in 
panel (a). They signal a first-order transition when m 
increases from zero to a finite value [see panel (a)]. The 
full (red) line is the numerical self-consistent result of 
Eq. (JTJ, and the dashed (blue) line is the approximate an- 
alytic result described above. The discrepancy between 
lines and symbols has a clear meaning. In order to obtain 
Eq. (JTJ we assumed that a second- order transition would 
take place. This is not the case, and the system under- 
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FIG. 3: (Color online) Hartree-Fock bands for j (full lines) 
and I (dashed lines) spin polarizations. 



goes a first-order transition for smaller U values. There 
are clearly two different regimes: one for Sn < 10~ 4 , 
where the dependence of Sn on U c is linear, and another 
for Sn > 10 -4 , where a plateau like behavior develops. 
This plateau has the same physical origin as the step like 
discontinuity we have seen in Fig. [T] (b). In the limit 
Sn — > we have not only [7 C — > 0, but also m — > and 
Am — > [see panels (a) and (b) of Fig. [2], implying a 
paramagnetic ground state for the undoped biased BLG. 

Figure [2] (d) shows also the effect of V on the T = 
phase diagram (the effect of t± being similar). Raising 
either V or t± leads to a decrease of the critical-?/ needed 
to establish the ferromagnetic phase for a given Sn. The 
order of the transition, however, remains first-order. We 
have observed that decreasing t± leads to a decrease in 
Am, and below some t± we can have Am < m. A sim- 
ilar effect has been seen when V is increased. It should 
be noted, however, that m and Am arc {/-dependent, 
meaning that, depending on V and t±, we can go from 
Am < m to Am > m just by increasing U . Irrespec- 
tive of V and t± we have always observed \5n\ < \An\: 
electron like carriers in one plane and hole like in the 
other. 

Understanding the asymmetry between planes. — The 
asymmetry between planes regarding both charge and 
spin polarization densities can be understood based on 
the Hartree-Fock bands shown in Fig. [3l Additionally, 
we note that in the biased BLG the weight of the wave 
functions in each layer for near-gap states is strongly de- 
pendent on their valence band or conduction band char- 
acter 0, [33l 34 1 . Valence band states near the gap are 
mostly localized on layer 2, due to the lower electrostatic 
potential —V/2. On the other hand, near-gap conduction 
band states have their highest amplitude on layer 1, due 
to the higher electrostatic potential +V/2. 

The case U < U c shown in Fig. [3] (left) stands for the 
paramagnetic phase. The values m = and Am = 
arc an immediate consequence of the degeneracy of | 
and I spin polarized bands. The presence of a finite 
gap, however, leads to the abovementioned asymmetry 
between near-gap valence and conduction states. As a 
consequence, a half-filled BLG would have n-2 = (4 + 
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FIG. 4: (Color online) Panels (a), (b), and (c) show the finite 
T solution for m, Am, and An, respectively, with T measured 
in K. Panel (d) shows the U vs T phase diagram. 



An)/ 2 e~/unit cell on layer 2 (electron like carriers) 
and n± = (4— An)/2 e~/unit cell on layer 1 (hole like 
carriers), with An 7^ 0. Even though the system is not at 
half- filling, as long as \8n\ < \An\ the carriers on layers 1 
and 2 will still be hole and electron like, respectively. 

Let us now consider the case U > U c shown in Fig. [3] 
(center). The degeneracy lifting of spin polarized bands 
gives rise to a finite magnetization, m =/= 0. Interest- 
ingly enough, the degeneracy lifting is only appreciable 
for conduction bands, as long as U is not much higher 
than U c . This explains why we have m w Am, as 
shown in panels (a) and (b) of Fig. [2]- as only conduc- 
tion bands are contributing to Am, the spin polarization 
density is almost completely localized in layer 1, where 
m\ = (m + Am)/ 2 ~ m, while the spin polarization in 
layer 2 is negligible, = (m — Am)/2 w 0. 

It is only when U 7$> U c that valence bands become 
non-degenerate, as seen in Fig. [3] (right). This implies 
that near-gap valence states with f and | spin polariza- 
tion have different amplitudes in layer 2. As the valence 
band for J, spin polarization has a lower energy the near- 
gap valence states with spin j have higher amplitude in 
layer 2 than their spin | counterparts. Consequently the 
magnetization in layer 2 is effectively opposite to that in 
layer 1, i.e., Am > m, as can be observed in panels (a) 
and (b) of Fig. H 

We note that the cases U >U C and U U c are param- 
eter dependent. The valence bands can show an apprecia- 
ble degeneracy lifting already for U > U c , especially for 
small values of the t± parameter. In this case the mag- 
netization of the two layers is no longer opposite, with 
Am < m. This can be understood as due to the fact that 
as t± is decreased the weight of near-gap wave functions 
becomes more evenly distributed between layers, leading 
not only to a decrease in Art but also in Am. 

Finite temperature. — Now we describe the phase dia- 
gram of the biased BLG in the T vs U plane. This is 
done in Fig. [5] for Sn = 5 x 10^ 5 e~/unit cell. For 
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T = — 1.1 K we studied the dependence of m, Am and 
An on the interaction U . First we note that the min- 
imum critical-?/ is not realized at T = 0. There is a 
reentrant behavior which is signaled by the smallest U c 
for T = 0.06 ±0.02 K. For temperatures above T « 0.1 K 
we have larger U c values for the larger temperatures, as 
can be seen in panel (a). The same is true for Am in 
panel (b). As in the case of Fig. [21 the value of Am, at 
a given T and U , is larger than to. Also the value of 
An, shown in panel (c), is larger than 5n. Therefore we 
have the two planes presenting opposite magnetization 
and the charge carriers being hole like in one graphcnc 
plane and electron like in the other. In panel (d) of Fig.Q] 
we present the phase diagram in the U vs T. Except at 
very low temperatures, there is a linear dependence of U c 
on T. It is clear that at low temperatures, T ~ 0.2 K, 
the value of U c is smaller than the estimated values of U 
for carbon compounds 

Disorder. — Crucial prerequisite in order to find ferro- 
magnetism is a high DOS at the Fermi energy. The pres- 
ence of disorder will certainly cause a smoothing of the 
singularity in the DOS and the band gap renormaliza- 
tion, and can even lead to the closing of the gap. We note, 
however, that for small values of the disorder strength the 
DOS still shows an enhanced behavior at the band gap 
edges 37]. The strong suppression of electrical noise in 
BLG [38( further suggests that in addition to a high crys- 
tal quality - leading to remarkably high mobilities [391 ] - 
an effective screening of random potentials is at work. 
Disorder should thus not be a limiting factor in the pre- 
dicted low density ferromagnetic state, as long as stan- 
dard high quality BLG samples are concerned. 

Let us also comment on the next-nearest intcrlayer- 
coupling 73 , which in the unbiased case breaks the spec- 
trum into four pockets for low densities [40j | . In the biased 
case, 73 still breaks the cylindrical symmetry, leading to 
the trigonal distortion of the bands, but the divergence 
in the density of states at the edges of the band gap is 
preserved [37] . Therefore, the addition of 73 to the model 
does not qualitatively change our result. 

Conclusion. — We have found that in the ferromagnetic 
phase the two layers in general have opposite magneti- 
zation and that the electronic density is hole like in one 
plane and electron like in the other. We have also found 
that at zero temperature, where the transition can be 
driven by doping, the phase transition between param- 
agnetic and ferromagnetic phases is first-order. 
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